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Abstract 

The upper-branch linear and nonlinear stability of compressible boundary layer flows 
is studied using the approach of Smith & Bodonyi (1982) for a similar incompressible 
problem. Both pressure gradient boundary layers and Blasius flow are considered with 
and without heat transfer, and the neutral eigenrelations incorporating compressibility 
effects are obtained explicitly. The compressible nonlinear viscous critical layer equations 
are derived and solved numerically and the results indicate some solutions with positive 
phase shift across the critical layer. Various limiting cases are investigated including the 
case of much larger disturbance amplitudes and this indicates the structure for the strongly 
nonlinear critical layer of the Benney-Bergeron (1969) type. Finally we show also how a 
match .with the inviscid neutral inflexional modes arising from the generalized inflexion 
point criterion, is achieved. 


‘Permanent address: -Mathematics Department, Exeter University, Exeter EX4 4QE, England. 
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§1 Introduction. 

The motivation for the present work arises from the need to consider the stability of 
compressible boundary layer flows but within the modern high Reynolds number asymp- 
totic framework. Many of the earlier classical approaches to the stability of incompressible 
boundary layer flows (see for example Reid (1965)) have now been largely superseded by the 
more self consistent asymptotic methods based on triple-deck or multi-deck ideas. Smith 
(1979a, b) was the first to apply such techniques to stability problems and he showed how 
the triple-deck structure governs the lower branch stability of boundary layer flows. This 
then allowed for weakly non-linear effects, strongly nonlinear effects, effects of the non- 
parallelism of the basic flow to be assessed systematically. Subsequent papers have applied 
the same basic concepts but with more complicated structures to consider the upper-branch 
stability of incompressible boundary layer flows, see Smith & Bodonyi (1980), Smith & 
Bodonyi (1982), Bodonyi, Smith & Gajjar (1983); the lower- and upper-branch stability 
of three-dimensional boundary layer flows, see Stewart Sz Smith (1987), Bassom & Gajjar 
(1988). 

As far as the stability of compressible boundary layer flows is concerned there have been 
several contributions by Lees & Lin (1946), Lees & Reshotko (1962), Reshotko (1962), 
albeit on classical lines. Pate & Schueler (1969) and Kendall (1975) present the relatively 
few experimental studies of compressible boundary layer transition. Most of this and sub- 
sequent work is extensively reviewed by Reshotko (1976), and Mack (1984,1986). However 
an important distinction between compressible and incompressible boundary layer flows 
is the existence of unstable linear inviscid modes satisfying the generalized inflexion point 
criterion, see Lees & Lin (1946), or Mack (1984). This would suggest that many of the 
linear results obtained by Lees & Lin are still valid. 

Our main objective in this paper is to extend the work of Smith & Bodonyi (1982) who 
considered the upper-branch stability of incompressible boundary layer flows, to compress- 
ible flows. We address upper-branch stability rather than lower-branch stability because 
the disturbance structure for upper-branch modes is more closely related to that for the 
inviscid modes than is the triple-deck structure. The importance of the upper-branch and 
inviscid modes is highlighted by the fact that these modes have shorter scales and higher 
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growth rates than the lower-branch modes. For a discussion of some aspects of the lower- 
branch linear instability properties of compressible boundary layers see the work by Smith 
(1987). Malik (1987) presents some recent computational results for hypersonic boundary 
layers. 

There have been very few studies of nonlinear effects in the stability of compressible 
boundary layer flows and in this paper we address at least some aspects of the nonlinear 
theory relevant for modes on the upper-branch. We show that the neutral modes are now 
governed by the properties of the compressible critical layer equation. The present work 
provides the foundation for studying nonlinear inviscid modes of either the generalized 
inflexion point type or the Benney-Bergeron (1969) type. In fact the scalings and structures 
for studying the latter can be deduced from the limiting cases that we consider. 

Starting with the Smith & Bodonyi (1982) paper, hereafter referred to SB, the most 
obvious strategy for studying the effects due to compressibility is to gradually increase 
the Mach number until significant changes from the incompressible theory occurs. For 
general pressure gradient boundary layers we find that the first major contributions due 
to compressibility occur when the Mach number becomes 0(1), for flows with and without 
heat transfer. In the outermost potential flow region the controlling equation is now the 
Prandtl-Glauert equation. Also if there is heat transfer then the density and temperature 
disturbances become singular in the vicinity of the critical layer. This then leads to an 
additional contribution to the disturbances shearing stress which balances the usual con- 
tribution at the wall and gives rise to the dispersion relations. The total contribution to 
the disturbance shearing stress at the critical layer is now proportional to the quantity 
D(pbDUb) occurring in the generalized inflexion point criterion. Here D = d/dY , where 
Y is the boundary layer coordinate, and Ub,Pb are the basic streamwise velocity and 
density respectively. For linear neutral modes to exist we require that the flow must be 
subsonic and that D(pbDUb) must be negative. Expressions for the neutral frequencies 
and wavenumbers are obtained explicitly for both the linear and nonlinear theories. In 
the nonlinear theory these neutral modes depend on the amplitude of the disturbance 
indirectly through the phase shift across the critical layer, as in SB, but now the phase 
shift has to be determined by solving a generalization of the Haberman (1972) equation 
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heat transfer. The nonlinear theory for pressure gradient boundary layers is considered in 
section 3 where we also derive the compressible critical layer equation. The properties of 
this equation are studied in section 4 and finally we conclude with some brief additional 
comments in section 5. 
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§2 Basic Equations and Linear Theory. 

The Navier-Stokes equations are nondimensionalized as follows with the suffix oo denot- 
ing local free stream values and an asterix denoting dimensional quantities. We introduce 
non-dimensional cartesian coordinates (z*,y*) = L(x,y), corresponding velocity compo- 
nents (u*,u») = Uoo(u,v), pressure p* = PooU^p, time t* = (L/U^t, density p* = p^p, 
temperature T!* = T^T, coefficient of viscosity p* = pooP, coefficient of bulk viscosity 
p[ = p'ooP 1 . This gives the non dimensional equations; 

(continuity:) 


d{pu) d(pv) 
dt dx dy 


= 0 , 


(2.1a) 


(momentum equations:) 


dufc 

dt 


+ 


Du/ 

dx. 


dp 1 
dxp Re 


d . , 2 d ‘ 

aT^ ~ 3' ,)v -“ + J 


(2.1b) 


(equation of state:) 


7 A&P = pT, 


(2.1c) 


(energy equation:) 



dp dp 


dp dp ' 

p 

dt Ua dx a 

= IP 

dt Ua dx a 


+P 


(7-1) 

Re 


1 | f l ^ \f ^ u o ^2 

-pe a pe a p + (p - 


+ 


ReM ^ &r a 



(2- Id) 


where 

du a dup 

tap = h ~ — • 

dxp dx a 

The nondimensional constants appearing in (2.1) are Re — UooL/v <*, the Reynolds 
number, a p the Prandtl number, 7 the ratio of specific heats, M 0 0 = Uoo/q the local Mach 
Number with q being the speed of sound. The characteristic lengthscale L can be taken 
to be the local station, say xo. We also assume that the Chapman viscosity law p = CT, 
where C is Chapman’s constant, holds. The Reynolds number Re is taken to be large 
throughout and we set the Prandtl number to be equal to unity. 
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§2.1 Pressure gradient boundary layer flows with heat transfer and insulated 
wall conditions. 

The basic boundary layer flow is given by, with y = Re~?Y, 

u = U b (x,Y), v = Re~W B (x,Y), p = p B (x,Y), 

T = T fl (x, Y) = 1/pb , P — Pb = 1/7 (2.2) 


For general pressure gradient boundary layers the following properties of the basic flow 
also hold, 


U B ,PB,T B -+ 1 as Y — ► 00 , 

U b ~\iY + \ 2 Y 2 + ... as Y-* 0+, 

Pb ~ Ro 4* R\Y + . . . as Y — ► 0+, 

Tb ~ So 4* SiY 4- . . . as Y — ► 04*) (2-3) 

where the coefficients Ai , A 2 , etc. depend on ®o and the Mach number. For insulated wall 
conditions Ri = Si = 0, and if there is heat transfer then Si ^ 0. The curvature term A 2 
depends on the pressure gradient and heat transfer. For further details of the basic flow 
see Stewartson (1964). 

We next introduce infinitesimal disturbances of size 8 say, (8 <C 1), and consider the 
stability of the basic flow. Since our interest is centered on the upper-branch stability of 
compressible boundary layer flows, the most obvious starting point is the work of Smith ic 
Bodonyi (1982). By gradually increasing the Mach number, we find that the first significant 
changes due to compressibility occur when M becomes 0(1). This is true also for Blasius 
flow with heat transfer ( with zero pressure gradient but ^ 0), although for Blasius flow 
with insulated wall conditions (Ri = \ 2 = 0) the significant changes due to compressibility 
arise much earlier. This special case is treated separately later in this section. 
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So taking the Mach number to be 0(1), it is appropriate to consider the scalings and 
disturbance structure for infinitesimal perturbations of the basic flow. The basic structure 
is as set out in SB, see also Fig. 1, and this continues to hold for compressible flows 
provided we consider wavenumbers of O(ile^) and frequencies of O(ReS), relevant to the 
upper-branch. There are five distinct zones Z1-Z5; the outer potential flow region Z1 of 
thickness 0(Re~&), the main part of the boundary layer Z2 , a thinner inviscid region Z3 
of thickness 0(Re~Ti) containing the critical layer; zone Z4 the critical layer of thickness 
0(Re~s $ ), and finally zone Z5 the wall layer of thickness 0(Re ~ & ). As in SB we consider 
disturbances proportional to E = exp[i(aX — or)] , so that 

(' u,v,p,p,T ) = (u B ,v B ,p B ,p B ,T B ) + 6{(u,v,p,p,9)E + c.c) + ... , (2.4) 

where c.c denotes the complex conjugate and x = e 5 X,t = e 4 r, and with the wavenumber 
a and frequency a expanded as 

a = a 0 + -f- . . . , a = a 0 + ta x + (2.5) 

Here e — Re~& is a small parameter. For neutral stability a 0 ,<xo are taken to be real. 

The major differences here from the work of SB arise because of the additional contri- 
butions due to the density and temperature. The disturbances now have the expansions 
in the various regions as: 

[euo + . . . , ei’o + ... , epo -j- ... , epo + • • . , + • • • ] 

in Z1 where y = e 5 y, 

[no + £U i + . . . , evo + e 2 iq + . . . , epo + e 2 p i + • • • , 

Po + e Pi + • • • , #o + + • • • ] 

in Z2 where y = e 6 Y, 

[i/°) + eu^ + . . . , + e 3 v^ + . . . , ep^ + e 2 p^ + . . . , 

p(°) + ep( x ) + . . . , 0(°) + efld) + . . . ] 
in Z3 where y = e 7 Y, 

[fi( 0) + . . . , e 3 t) (0) + . . . , ep< 0) + . . . , epW + . . . , e^ 0) + . . . ] 
in Z5 where y = e 8 y, 
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together with appropriate expansions for the basic flow using (2.2),(2.3). 

In the outermost potential flow region the leading order equation is now the Prandtl- 
Glauert equation, 

“ a o(! ~ M lo)Po = °- 

Anticipating the requirement that M < 1 for neutral solutions and applying a bounded- 
ness condition at y = oo gives 

Po = Poe~ aaU>i> , vo = -iPoue ~ aoU,f> , u 0 = -po, 

7 Po = Po , Oo = yM^po - po, (2.7) 


with u = (1 — M £,)^. Also P 0 is an unknown constant. 

Solutions in the other regions follow in a straightforward fashion and the important ones 
are given by 


u 0 = Ao 


dU B 
dY ’ 


v 0 = -iAoaoUs, 


- A dpB 

P 0 AO Qy } 


Oo = -po— , Po=Po, u (0) = A 0 A!, 
PB 


u(°) = 


l ^^-iaoAo{\iY-co), p (0) = 


R lV W 


P 0 Ai 


icto^F - c 0 )’ 


0(°) = 


SopW 

Po ’ 


p < 0) = a, 


„(■) _ fi (>) _ iaoj4 ° c ° 

A 




( 2 . 8 ) 


Here we have put the wavespeed Co = cq/qo, Aq is an unknown constant, and £ = (AjF — 
co). Also contains terms which are regular as ( -* 0. From (2.8) it is seen that as 
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the critical level is approached (£ — » 0+), the leading order solutions for the density and 
temperature fluctuations in zone Z3 become singular- if we have heat transfer, i.e. if Ri ^ 0. 
See also Lees 8z Reshotko (1962), and Mack (1984). For compressible flow the coefficient 
of the logarithmic term in the solution for t/ 1 ) now has a contribution proportional to 
( Ri/Ro ) arising from the variation of the basic density or temperature. It can be seen 
that the total contribution is now proportional to (pb^y~) |y=o- 

Matching between zone Z1 and Z2 and requiring that the normal velocity fluctuation 
goes to zero at the wall then gives from (2. 7), (2. 8), 


and 


Po = a 0 Ao/u> 


Pq — j4o c o-Rq^i. 


Hence we obtain the relation connecting the wavenumber and wavespeed 


£*o = cqRo-^iw = coRoA^l — M^,) 2 . 


(2.9) 


If M > 1 then with the solution for outgoing waves from the Prandtl-Glauert equa- 
tion, the details in the other regions Z2-Z5 remain largely unaffected. However when the 
matching is performed, because of an extra i factor now introduced in (2.9), this shows 
that no neutral solutions for > 1, with these scalings, exist. With = 0 in (2.9) 
this reduces to the relation obtained by SB. 

The singularity in the density, temperature and velocity fluctuations in zone Z3 requires 
a thinner region the critical layer of thickness 0(e 3 3). However for conciseness we summa- 
rize the main result that the solutions in Z3 continue to hold for £ < 0 provided we take 
ln£ — > (In |£| — in) for £ < 0. The details of the linear critical layer may be deduced from 
the nonlinear theory in section 3 below as a limiting case for small amplitudes. The wall 
layer is essentially the same as that in SB and gives the condition 
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where now 


m = 


Ro^o 

Hb 


i 

i 


and hb = n(So) is the coefficient of viscosity at the wall. 
Matching with Z3 therefore yields the second relation 


OfOCQ [2A2 Ri 


<*qCqAi ( hb \ 


( 2 . 10 ) 


Ai [Ai Ro\ K ’ <x 0 \ 2 Ro<to) ' 

Combining this with (2.9) we finally obtain the linear neutral eigenrelations for the upper- 
branch stability of compressible boundary layer flows, 


and 



(2.11a) 


(2.11b) 


(2.11c) 


For the above relations to hold we require, in addition to the condition for subsonic flow, 

2A R 

that the quantity x = + -5-) be negative. 

Ai Ro 

If = 0, and we set fiB = Ro = 1,-Ri =0, then the results of SB are obtained. (The 
latter conditions are necessary for consistency with the work of SB where a uniform basic 
density and temperature field is taken.) We note that the (1 — M^,)» factors cannot be 
rescaled out of our problem because R\ depends on 

Some interesting limiting cases of the relations (2.11) can be examined further. Firstly if 
X — > 0— with the Mach number fixed, then from (2.11ab) c 0 ~ (— x) - ® an d no ~ ( — x) - ^- 
Hence formally if x becomes 0(e 3 = Re~*) then since the effective wavenumber a is 
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a = e 5 a 0 + . . . , we obtain 

eco ~ 0(1), a ~0(i?e^) and a ~ O(e~ 4 cro) ~ O(Re^). 

This implies that the critical layer moves away from the wall and into the main part 
of the boundary layer. The limiting neutral structure is then precisely that of classical 
form with inviscid scalings, i.e. the streamwise and lateral scalings are comparable to the 
boundary layer thickness. Also the condition x 0— implies that the above neutral modes 
match on to the neutral inflectional modes which come from the generalized inflexion point 
criterion. The above situation is similar to that occurring in the upper-branch stability 
of three-dimensional incompressible boundary layer flows, Bassom & Gajjar (1988). It 
should be noted however, (as a referee has kindly pointed out ), that by varying x the 
properties of the basic flow are being changed. In Bassom & Gajjar the Rayleigh scalings 
were obtained by changing the orientation of the three-dimensional disturbance. 

As an example if we consider Blasius flow with heat transfer, then 



where D is some positive constant and Su, is the ratio of the prescribed wall temperature to 
the recovery temperature, Stewartson (1964). The condition that x < 0 means therefore 
that the neutral modes exist only for subsonic flow with cooled walls (s«, < 1). As we let 
s w — ► 1— then x ~ ^ 0 — ? an d the neutral structure shrinks and we go over to the inviscid 
scalings. If we include small pressure gradients then there are some values of s w (s w < s Wc 
and s Wc >1) for which the neutral solutions exist for heated walls also. Again in the limit 
s w — » s Wc ~ we obtain the inviscid scalings. 

The second limiting case concerns the limit M — ► 1 — . Prom (2.1 lab) we find 

c 0 ~(l — M£,)-A, a 0 ~(l — M^o)£. (2.12) 

Formally if (1 — M becomes 0(e 12 = i?e -1 ) then again the critical layer moves away 
from the wall. But now the streamwise lengthscale increases and the frequency decreases 
to both become 0(1). Since also from the solutions in the outermost region Z1 e 5 y ~ 
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e 5 a^' 1 (l — M^) - * becomes O(Re^), then this suggests a three tiered structure with a 
streamwise lengthscale of 0(1), and regions of lateral extents 0(Re~i), 0(1) and 0(Re*). 
The implications and importance of this transonic regime require further study. Professor 
F.T. Smith (private communication) has indicated that a similar situation arises in the 
transonic regime for the lower-branch modes, see also Smith (1987). 

Lastly we note that if we increase the disturbance size 6 and allow for nonlinear effects 
then, with 0(1) Mach numbers, and as in SB the critical layer becomes fully nonlinear 
when 8 becomes 0(e%). The structure and analysis for the nonlinear compressible critical 
layer is considered in section 3. 

Before that however we return to the case of Blasius flow with insulated wall conditions. 
It is clear from (2.11) that with A 2 = Ri = 0 the above analysis breaks down. The main 
reason is that the effects due to compressibility start coming in for much smaller Mach 
numbers, because the generalized inflexion point is in the main part of the boundary layer 
for 0(1) Mach numbers but retreats to the wall as — ► 0. 

§2.2 Blasius flow with insulated wall conditions. 

When M becomes 0(Re~^) then the neutral eigenrelations for the upper-branch of 
Blasius flow with insulated wall conditions will be significantly modified from the incom- 
pressible case. A rough argument to show this is as follows. If we consider the neutral 
modes satisfying the generalized point criterion, for Blasius flow with insulated wall con- 
ditions, then these satisfy the condition 

W i PB lsf) = °' Ub = C ’ at Y = Y °' ( 2 -l 3 ) 

where c is the phase speed of the wave, and pg, Ug are the basic density and streamwise 
component of velocity. 

The basic velocity and density are given by 

Ub = fin), T b = 1/pb = 1 + A&U - (2.14a) 
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where f(rj) is the Blasius function satisfying 


/"' + //" = o, f(0) = /'( 0) = 0, /'( oo) = 1, 
and fj is related to Y by the Dorodonitsyn Howarth transformation 


(2Cx 0 )?fj = [ 
Jo 


p B dY. 


For small tj, (2.13) yields 


where 


Hence 


-rj [2(7 - + 12A 4 ij] ~ 0, 


f ~ A^ + A 4 j? 4 + . . . as r] — » 0. 


rj ~ 0 or r) ~ 


(7 - lyxiMi 
6A 4 


(2.14b) 


(2.14c) 


(2.15a) 


(2.15b) 


Now for small c the critical layer position is given by Y ~ c/ Aj , and since on the upper- 
branch of incompressible Blasius flow c ~ 0(Re ^), Smith & Bodonyi (1980), this sug- 
gests looking at the scaling 0(Re~&). 

If we set €i = Re~?5 , and , then the basic flow has the following properties, 


Ub(x,Y) ~ Ubo + £iUbi + • • • 

Vb = 0(e}°), PB — 1 + + • • • > 

Tb — 1 + ^ 1^51 + • • • , PB = e 1 1 pB = 3 v 2 1 

7^00 

and close to the wall with y = ej 1 V', 


Ub ~ t\ a \\Y + t\a.2\Y + Cia 3 iF -f- e{(a 4 iF + a 4 3F 3 + a 44 F 4 ) + . . . , 

PB ~ 1 + €i 6io + €5620 + ^1(^30 + & 32 F 2 ) + . . . , 
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where the constants a{j,bij can be obtained explicitly from (2. 14), (2. 15), and in particular 


an = Aj = \ 1 /(2Cx 0 )$, 6 a 43 = Aj( 7 - 1 

a 44 = A 4 = —Xj/24(2Cx 0 ) 2 , 2 b 32 = A*( 7 - 1)M^. (2.16) 

The linear disturbance structure for incompressible flows is set out in Smith & Bodonyi 
(1980), see also Fig.2, and this work can be extended to compressible flows in a straight- 
forward fashion, using the properties of the basic flow as given above. For brevity however 
only the main points are outlined. The wavenumber and frequency now expand as 

9 9 8 

a = Re™ a ~ Re™do + Re™at\ + . . . , 

8 8 

& = Re 20 a ~ Re 30 ( 7o + ... , 

respectively, and Co = (Tq/oq is the phase speed. We set x = efXi and t = t\r\ and look 
for disturbances proportional to E\ = e , ^ aXl ~ erTl K Then in the main part of the boundary 
layer, the disturbances have the expansions 


u — u 0 + t\U\ + ... 

v = eiv 0 + €ji>i + . . . 

p = eipo + 4pi + ..., 

P — € iPo + e iPi + • . . , 

9 = eiOo + + . . . , 

Substitution into the Navier-Stokes equations produces the leading order solutions 


i / 0 = AqUboy , 


Po — —9 q 


Vq = — iOfo^o^BOj 

dpBi 


Po = -P 0 , 


■Aq- 


dY 


where Aq,Pq are constants. By writing down similar expansions for the outermost potential 
flow region and matching with the solutions above gives the relation, 
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P 0 — a oAo- 


as in the incompressible case. 

Next in the thinner inviscid region which is part of the main boundary layer containing 
the critical layer, with y = the disturbances expand as 


u = u(°) -f + . . . , 

v = e?u<°> + + e{v< 3 > + . . . , 

P = e iP (0) + e?P (1) + ...» 

P = e ?P (0) + «?P (1) + • • • , 
e = ^( 0 ) + e 3 6 (l) + 

The leading order solutions for td°\ v (0 \p (0 \ are as in SBl and produce the relation 


m 0 P<°> 

an 


= i<7o.Ao> 


after requiring that the normal velocity fluctuation tends to zero at the wall. The solution 
for the density disturbance yields, 


P (0) ( Aoa u 2b 32 Y 2 
7 PB ( a u Y-c Q )' 


Again as for the previous case studied, the density and temperature fluctuations become 
quite large in the vicinity of the critical layer. The higher order solutions can be obtained 
by successively solving the equations at each order. As in SBl it is necessary to proceed 
to O(ef) in the expansions for v, where we find that 

u ( 3) = PLT - 4-(« n-Fo + c 0 Pi)(a„f - co)ln(a„y - c„), 
a 11 
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where PLT contains algebraic terms which are regular as Y — ► co/au, and 


Fo = 2iAoaoaii(a n b 3 2 + 3043 ), Fi — 12mo a n a 44-^o- 

The Fi contribution to the logarithmic term is as in the incompressible case but the 
additional contribution due to the Fo term, arises from compressibility effects. 

The critical layer needs to be introduced to smooth out the singularities in the density, 
temperature and velocity fluctuations, but the ensuing result is the familiar *7r’ jump in 
the logarithmic term across the critical layer. The leading order wall layer solutions are as 
in SB1. After completing the matching, we finally obtain the eigenrelations, 


a 0 = coan, 


(2.17) 


and 


2 cl 


~(— fl")(& 32 a n + 3an<l43 + 600044) = 


a 0 


*11 


{2alcl)h 


Substituting for 632,043 and a 44 , from (2.16), in terms of the scaled Mach number gives, 

(2.18) 


H|=[At(7-l)*i+6c„A 1 ]=--|2 J 

A i (2<*o c o 


) 5 


Hence for the neutral modes to exist we require also from (2.18) that 


6c 0 > 


a {(7 - mi 


The eigenrelations are plotted in Fig. 3a. If = 0 then (2. 17), (2. 18) reduce to the 
expressions obtained by Smith & Bodonyi (1980). Clearly the neutral solutions exist only 
for a certain range of wavenumbers and frequencies. If becomes large, then from 
(2.18) we get 


co 


6A. 




(2.19) 


where 0 < cj << 1. Hence for increasing Mach number the critical layer moves away from 
the wall. When M becomes O(Re^o) formally, i.e. the unsealed Mach number becomes 
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0(1), and from (2.17) and (2.19) it is seen that the wavespeed becomes 0(1), and the 
wavenumber and frequency become O(Re^). So again the inviscid scalings are recovered. 
Prom (2.15a), (2.19) it is clear that in this limit a match with the neutral inflexional mode 
1 for low 0(1) Mach numbers is achieved. If we replace cq by c and M by Moo in (2.19) for 
0(1) Mach numbers, then the predicted asymptote agrees well with the computed phase 
shift c from the generalized inflexion point criterion, see Fig. 3b and also Fig. 9.2 in Mack 
(1984). 
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§3 Nonlinear compressible critical layer equations. 

As in SB when the disturbance size S becomes 0(e 5 ) then nonlinear effects in the 
critical layer become important and cannot be neglected. For Blasius flow with insulated 
wall conditions the scalings and disturbance size to produce a nonlinear critical layer are 
somewhat different and are not considered here, but see Cole (1989) for a discussion of 
this case. 

Since we are interested in wavenumbers of 0(Re&) and frequencies of 0(Res ) we set 

d . d , d 

— ,)gj + — 

d ^ 

where as in SB the slower d/dx variations do not play any significant role, and can be 
ignored. The constants a-o, oi , 00, 01 ... , are taken to be real, so that we are considering 
neutral disturbances and the wavespeed Co = 0o/c*o- 

Before proceeding with the critical layer analysis, some brief details of the solutions in 
the other regions are necessary. 

In zone Z2, the main part of the boundary layer, with y = e 6 Y, the total expansions for 
the flow quantities now take the form, 

u = Ub(x, Y) + e*u Q -f e^ui + . . . , 

10 13 

v = + € 3 Vi + . . . , 

P = PB + + e^Pi + • • • , 

p = p B (x, Y) + e3/? 0 + pi + • • • , 

T = T b (x, Y) + eh 0 + e 1 ?6i+.... (3.1) 


Substitution into the Navier-Stokes equations then produces the leading order solutions, 


- my 

u 0 AyJC) qy 1 


dA(X) TT ~ dp B 

Vq = -CXQ— Qjr— U B, Po = A ~gy~> 


a _ TbiQpb 

VO — A _ nv r > 

PB oY 


Po = Po(X). 


(3.2) 
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The unknown disturbance function A(X ) will be set equal to Ao cos X later, but for con- 
venience we assume the more general form for now. 

In the inviscid region Z3 containing the critical layer, with y = e 7 Y, the following 
expansions are implied from (3.1), (3. 2): 

u = e\i Y + e 2 A 2 Y 2 + + €*umf + e 3 A 3 F 3 + + . . . , 

P = PB + ^ P (0) + e^p (1) + . . . , 

p = R 0 + eRiY + 6 *pmf + e 2 R 2 Y 2 + esp (0) + . . . , 

T = S 0 + eS\Y + eh MF + e 2 S 2 Y 2 + eM°> + .... (3.3) 


The suffix MF denotes mean flow terms (independent of X) which arise because of the 
interaction of the nonlinear critical layer with the outer flow. 

Substitution into the Navier-Stokes equation produces the solutions 


u (0) = iAj, u (0) = -a 0 Ax*iY, p (0) = P (0) , 


where the relation 

a 0 3p (0) - 

P 0 Ai dX ~ aoAxCo ’ 

is required in order that the solution tends to zero at the wall. The density disturbance 
p ( ' ’ J becomes singular at the critical layer and this gives an additional contribution to the 
logarithmic term at higher order. In fact we have 


(i) - _ 


ir ’ = 


Ro Ai 


dpW dp ^ du^ 

«o o — + + Po(coai — cTi)- 


dX 


dX 


dX 


-Fo(n+^(HO-A) 

Ai Ai 

Ai Ax 

where £ = (AiF — c 0 ), Ai is an unknown function of X and 


(3.5) 


F i = 


QqA 2 Ax 
Ai 


Fi = 


agRi cq Ax 
Ro 
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The solution for can then be obtained from the continuity equation 


du (1) du (0) du (1) n 

a °~dx JrC ‘'-dx + ~W- a - 


f 23 

Hence in the critical layer Z4, where y = + e~* 1 ,, we have from (3. 3), (3. 4) and (3.5), 

:i«(D 


e2A 2 C 0 , A.-.( 1) , ,#f,(2) , M (3) 


u = ec 0 + e*\ir) -\ ^ 2 — ~ + e s £r > + e 5 tr ’ + e — -ip + e 3 u ( } + — 

A x Aj 

v = e*®<o> + ex 5 0) + «¥«<*> + e %v c (X) + e^ 3 > + . . . , 

P = Ro + + e5/5 (0) + e*p {1) + . . . , 

Ai 

T = S 0 + + e!0<°> + . . . , 

p — p B + p^ + p^ + e 1 * p ^ p^ + 


(3.6) 


Substitution into the Navier-Stokes equations then gives, 

= — a 0 \iAxV, = AAi , 

2 c 0 dv (2) _ 
dr, ~ A: ’ dr, ~ 

At the next order, the governing equations are 


v<°> 


<*o 


dpW 


AjiZo dX 


du (2) 


Ro 


d , ./ 1U dv (3 * 

(a 0 u {3) +ai« (1) ) + 


dx 


dr. 


+ 


Hi CO 
•^1 


a 0 - 


dv 


dX 


+ 


dr, J 


x dp (0) ~ (0 )dp {0) 

+ A I ,a°-^ + »< ) -^- 


= 0 , 


together with 


Ho 


. du^ du^\ 

Al ’'( a “-sF + “ , lF ) + " o( “ 9F 


+ u< 2 > 


du^ 


dX 


7^)+ 




dr, 


dr. 


+ v™ 


dii^ 

dr, 


+ w (3) Ai 


+ 


Hi c 0 rrtl , ~( 0 ) ~(o) \ d P {3) , d , 5 “ (3 \ 

^-[0] + '’ ” Al = -““W + 
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and 


and finally, 


7 PB 


a 0 \iV 


dpW 

dX 




dpW 

di) 


fi B Ro <9 2 0 (o) 

ML dr l 


P {0) S 0 + R 0 fc 0) = 0. 


The pressure terms in (3.6) are independent of 77 to the order indicated. It is at this order 
that the first non-simple solutions are generated which smooth out the irregularities in p (0) 
and i>U), Using the previous solutions we obtain after some manipulation the nonlinear 
compressible critical layer equations, 


, „ a ,aa<»\ . „ „, 0) 8 2 fi< 3 > a 3 « (3) >2 dp<°> , , 

ooA.^Bo— (-^-) + B 0 » ( < 3 - 7a) 


where p^ satisfies 


x dpW ~ (0 )dp (0) 

<x 0 \m~r + v (0) -!— 


dX 


drj 


q av 01 n 

- W 5„-g^- = 0. 


(3.7b) 


We note that for a more general viscosity law, the Chapman constant C in (3.7a) is replaced 

^So) 

dT 

From (3.4)-(3.6) the boundary conditions are 


R\ Ac§ 

as 77—7 ±00, 

(3.8a) 

Ai 77 

,2A 2 , 

Aj + 

— as 77—7 ±00. 

Ro v 

(3.8b) 


^ ~ 2M n + H ± 

Also matching with Z3 gives 

= — (XqAxCq, 

so that if we set A = Aq cos X and introduce the normalizations 

77 = gZ, g 2 = ° ° , A c = pb/(9 3 XiR 0 (Xo), 


'1 


du^ c 0 A 0 ,2A 2 ~(o) n rr 

or) g A : i? 0 


(3.9) 
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then (3. 7), (3. 8) become 


and 


aa . v an . d 2 si Tff7 du . a 2 n. 

Z dX + Sm A dZ Xc dZ 2 ~ J ^ Z dX Ac dZ 2 


. an . a 2 n 0 
z dx +smX dz Ac az2 _0 ’ 


(3.10a) 


(3.10b) 


with the boundary conditions 


Cl~(l-J)Z+ < ^- + H ± 


as Z — ► ioo, 


„ „ cos 7C A + 

n ~ z + — =— + a* 


as Z — ► ± 00 . 


(3.10c) 


(3.10d) 


Here the parameter J = R\\\/{2\2Rq + AiiZi). The set (3.10a-d) constitutes an exten- 
sion to compressible flow of the Haberman (1972) equation. The nonlinear compressible 
critical layer equation differs from the usual Haberman equation because of the extra forc- 
ing term on the right hand side of (3.10a). The parameter J depends on Ri, the basic 
density gradient at the wall and A 2 the curvature of the basic flow. For problems involving 
heat transfer J is nonzero. Also as in the linear theory the above analysis applies equally 
to Blasius flow with heat transfer when J = 

The important quantity as far as the eigenvalue problem is concerned is the velocity 
jump and phase shift <f>, resulting from the set (3.10). 

If we put Q = dU /dZ, then from (3.10c), 

2 2 

U ~ (1 - J)— + H±Z + cos X In \Z\ + U±(X) as Z -* ± 00 . 

The phase shift <f> is then given by 

<f> = — [ sin x(U + — U~)dx. 

n Jo 

and as in Haberman (1972), SB, the following identity can be derived, 

<j> + J4 > h =2X c (H + -H~). (3.11) 

Here <^//(A c ) is the phase shift from the Haberman problem, with 7 = 0 above. 
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The problem (3.10) has to be solved numerically to determine <f>, but as far as the 
eigenvalue problem is concerned we note firstly that ln£ — ► (In |£| + i(f > ) for f < 0, in (3.5). 
Then since the leading order solutions outside the critical layer are determined from linear 
equations, the eigenrelations can be obtained as in section 2 with now — ir being replaced 
by <f>. 

This then yields the relations for neutral stability, 


co = 


4 


(1-JW*)A 


'2A 2 


+ 


RcJ 9 l 


l* G*s) 


<*0 = 


Ro\hl-Ml)k (» B \‘ 

\ 2 Rl) 


and 


<?o 


[Y2A 2 

Ri' 

\ ' 

[(a, + 

Roy 


i2oA?(l 

-Ml)* 

[ 72 A 2 

R,) 


lUr + 

RoJ 

<f> 

- 


(3.12a) 


(3.12b) 


(3.12c) 


i \2RlJ ' 

Clearly for neutral stability we require the condition M < 1, i.e. subsonic flow and 


O \ Tf 

<f>x = + £-) > o. (3.i3) 

Ai /l0 

Finally since the equation for the density disturbance II is the Haberman equation, 
(A + — A - ) is non zero, and hence there is a jump in the mean density across the critical 
layer. As the parameter A c goes to infinity we recover the linear results since <f> — > — 7T. 
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§4 Results and Discussion. 

The problem (3.10) was solved numerically using a Fourier series for II and Q, and solving 
a truncated system using second order central differencing. The discrete equations were 
solved iteratively using relaxation especially for the small A c values. The usual checks were 
made on grid size independence, the number of Fourier terms used etc. The phase shift 
<j> and the vorticity jump (H + — H~) are shown in Fig. (4ab) for several values of J. For 
J — 0 our results agree with those of Haberman (1972) and SB. Solutions were obtained 
for both positive and negative values of J. For J becoming large and negative, (j) becomes 
much more negative than the J = 0 case. For J positive however, for a fixed value of 
A c the phase shift is larger than the equivalent Haberman value. Our calculations also 
indicate that there exists critical values of J = J c > 0 such that for J > J c , the phase shift 
becomes positive for some values of A c . 

Clearly J < 0 implies either Ri/R 0 > 0 and x < 0, or R\/Rq < 0 and x > 0. Our 
results indicate that no solutions are possible with the latter case since from (3.13) we 
also require <f> > 0. The former case implies a positive basic density gradient (and hence a 
negative temperature gradient) at the wall, and applies only for accelerated flows. 

If on the other hand J > 0, then we require that R\ /Rq > 0 and x > 0 and hence <f> > 0, 
or Ri/Rq < 0 and x < 0 ( and </> < 0). The numerical results show that there are solutions 
with positive <f> if J > J c > 0 and 0 < A c < A cr u(J), where A cr a is the value of A c such 
that <f>(J, A c ) = 0. Since R\/Rq > 0 this implies heated walls only. The other case with 
<t> < 0 has solutions only for J < J c , \ cr u < A c , and requires cooled walls. 

The nonlinear neutral curves will differ considerably depending on the above cases. Some 
sample curves of the nonlinear neutral amplitude against frequency are sketched in Fig. 5. 
It can be seen in Fig. 5a that if J < J c then the behaviour of Ao is qualitatively similar 
to that for the incompressible (SB) case. If however J c < J and <j> < 0 then (Fig. 5b) the 
amplitude reaches a peak and decreases to zero as the neutral frequency increases. Finally 
if J c < J and <f> > 0, then Fig. 5c indicates the possible behaviour. We note that these 
latter solutions are not present in the linear case as the amplitude approaches zero. There 
is a critical value of the neutral frequency below which the nonlinear solutions disappear. 
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In the limit J — » oo, A c = 0(1), a limiting solution of (3.10) can be found in the form 


Cl = JCl\ + . . . , <f> — J(f>x + . . . , 


where Cl\ satisfies 

ZQ,\x + simftiz — A c f)i zz — ZYlx — ^ c Clzz (4.1a) 

with the boundary conditions, 

fti ~ —Z + H ± as Z — ♦ ioo. (4-lb) 

The phase shift from the solution of (4.1) is shown in Fig. 6 and this is in good agreement 
with the numerical solutions of the full equations, see Fig. 7. For |J| large the above 
predicts large mean vorticity gradients and phase shifts across the critical layer. 

The limit A c — > oo, | J\ — » oo, J = 0(A C ), can be analyzed in the same manner as the 
linear limit with J = 0(1). This then reproduces the linear result that <j> — ► — it as A c — ► oo, 
with J = 0(A C ). 

The last case which merits further consideration is the limit A c — > 0 when the critical 
layer becomes strongly nonlinear. Our numerical results indicate and also the correspond- 
ing incompressible limit suggests expansions of the form, 

ci = fio + Acfij + . . . , n = n 0 + A c rii + . . . . (4.2) 




Here No is a constant, and (4.4) expresses the fact that the density is constant inside the 
closed streamlines 77 , < 1, and is continuous across the boundaries, see SB for further 
details. 

After substituting (4.2) into (3.10) the equation for fio produces 

= JK'McosX + GoM, (4.5) 


where Go is an unknown function of 77 , but again determined from a viscous secularity 
condition at higher order. The equation for ill is 

5Qi r . an j d /rr dCl 0 

= -J sin A— — (Z — — ), 


dX 


drj* dij t dr)* 


with 


f X d 

ill = ~ — (2 *( 77 , - cost)>K')dt + K\{r)*) 77 * > 1, 

Jo & 1 ]* 


= Ni r)* < 1, 

where Ni is a constant and Ki is an unknown function of 77 ,. Hence using the above 
expression for III , and applying a periodicity condition on gives for 77 , > 1 , 

r2n I-Zi q2 x j 

— J / sin# / — —(2 7 ( 77 * — cosxi)? K')dxi dx 

Jo Jo dril 

r2ir 


+ - — / 25 ( 77 , - cosx)* — — dx = 0. 

5?7, Jo dll* 


After using the boundary conditions as 77 , —7 00 this now gives 


G'o = K' 


(1-7) + J- 


f2n COS X 

L r dx 

2(77, — COS X ) 3 

fo*( 77 * — cosx)5 dx 


(4.6) 


For 77 , < 1 the approach adopted by SB, Brown-Stewartson (1978), and Goldstein & 
Hultgren (1988) gives G' 0 = 0 inside the closed streamlines. From (4.6) G' 0 becomes 
logarithmically singular at the boundaries of the closed streamlines as 77 * —7 1 +, and we 
find that 

+ a, ,.-1 + .. (4.7) 
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Viscous layers are thus necessary near the boundaries of the cat’s eye to resolve the 
discontinuity in the vorticity gradient. Although the vorticity gradient is discontinuous 
also near the boundary 77* = 1 in the incompressible case, here there is a much stronger 
discontinuity because of heat transfer. Details of the viscous layer calculation are omitted 
here, but in summary we find that the phase shift from the viscous layer alone is 


In 16 1 

— A c 7r[7(-yg- + -) 




across each layer. The phase shift across the inviscid region, from 77* = 1 to 77* = 00, above 
the cat’s eye is 


- 27) - *(1 - y) + c< 2 > J ~Y J + !n(16)]. 


The same phase shift occurs below the cat’s eye. Here C ^ is the constant appearing in 
the incompressible theory, see SB, and 


c ( 2 >= ~ rv r 

2?r Ji J 0 


cos X 


(t?* - COS .T ) 2 


- dx drj *. 


Our calculations give 


C (1) = -5.516, C (2) = .194 . 


Thus the total phase shift is found to be given by 


<t>~ A c [C (1) (l -2J) + 2C (2 ) J-2ttJ] as A c -> 0. (4.8) 


The above calculation neglects the smaller contributions from the edges of the cat’s eyes 
at X = 0, 2ir. 

For J = 0 the above asymptote agrees with that in SB. For J ^ 0, the last term in (4.8) 
represents a contribution from the viscous layers which does not cancel out, in contrast to 
the incompressible case. 

Comparisons of (4.8) with the full numerical calculations are encouraging and the general 
trend for different values for J is correctly predicted, see Fig. 8. 


28 



For fixed Mach number and as A c — > 0 (3. 12), (4. 8) shows that the wavespeed and 
wavenumber increase and we recover exactly the scalings of SB and Bodonyi, Smith Sz 
Gajjar (1983) for the strongly nonlinear critical layer problem. 

Next for fixed properties of the basic flow, the numerical results show that there are 
cases when 4> — ► ±0 as A c — ► ±A cr j t (J), for 0(1) values of A cri( . Then (3.9) together with 

(3.12) shows that Ao ~ o 0 3 c^" 1 and c 0 ~ |<^| - 3,a 0 ~ Thus formally in the limit 

when Co becomes 0(e -1 ) the disturbance amplitude (eSAo) falls to become 0(Re~ 5). In 
this regime the Rayleigh scalings are obtained but with a viscous nonlinear critical layer 
and with zero phase shift across the critical for the leading order stability problem. 

As x ~ * 0, for A c = 0(1), the situation is different from the linear case. Since x —■ * 0 
implies |J| — + oo for Ri/Rq ^ 0, we have that \ Rifi/Ro to leading order. From 

(3.13) we need R\/Rq positive and so this limit applies for heated walls only. In order to 
retreive the Rayleigh scalings of the linear theory the limiting process A c — ► oo, J — > oo, 
but J — 0(A C ) is necessary. 

Finally the last case concerns the limit M — * 1 — . Again for fixed A c = 0(1), the 
wavespeed and the wavenumber become 0(1) when (1 — M '£,) ~ 0(l/i2e). But from (3.9) 
we have A ~ o 0 5 e 0 , so that in the above limit e*A becomes 0(1). Thus in the transonic 
regime, for the nonlinear problem, the disturbance amplitude becomes 0(1) comparable 
with the basic flow. The earlier comments regarding the three-tiered structure still hold 
and again this special case deserves further study. 
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§5 Summary. 

We have shown how the work of SB,SB1 can be extended to compressible flows. In 
particular the results of SB continue to hold for compressible flows with insulated wall 
conditions until the Mach number becomes 0(1), but for Blasius flow with insulated wall 
conditions the effects of compressibility start coming in when the Mach number reaches 
an 0(Re~ w) value. In the compressible theory the density and temperature fluctuations 
play an important role in determining the properties of the nonlinear compressible critical 
layer and hence the neutral modes. 

Several results with possible important practical implications are that the density and 
temperature disturbances are quite large near the critical layer, and in the nonlinear theory 
the interaction of the critical layer with the outer flow produces mean flow distortions in 
both the temperature and density fields. Also with increased disturbance amplitudes the 
nonlinear critical layer takes on a form similar to that of incompressible flows with open 
and closed streamlines. 

We have seen also how the classical inviscid scalings for the neutral inflexional waves are 
retrieved as the quantity D(pbDUb) goes to zero. The nonlinear disturbance structure for 
this limit, and also that for the strongly nonlinear critical layer, is currently being inves- 
tigated. This should provide further new results on the nonlinear aspects of compressible 
boundary layer stability. In the strongly nonlinear case, as in SB and Bodonyi, Smith & 
Gajjar (1983), the influence of the transverse pressure gradients across the critical layer is 
not negligible, and in addition the effects of the varying density and temperature also has 
to be taken into consideration. How this affects the nonlinear neutral states remains to be 
seen. 

Finally we emphasize that although equilibrium critical layers are a special case, the 
current work provides the scalings and structures for investigating the possibly more im- 
portant problem of spatially and temporally growing modes. 
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FIGURE 1. - SKETCH OF THE VARIOUS REGIONS FOR GENERAL BOUNDARY LAYER 
FLOWS. Z1 IS THE OUTER POTENTIAL FLOW REGION. 12 THE MAIN PART OF 
OF THE BOUNDARY LAYER. Z3 THE THINNER ADJUSTMENT REGION CONTAINING 
THE CRITICAL LAYER ZA, AND Z5 IS THE WALL LAYER. 



FIGURE 2. - SKETCH OF THE DIFFERENT REGIONS FOR BLASIUS FLOW WITH 
INSULATED WALL CONDITIONS. Z1 IS THE OUTER POTENTIAL FLOW REGION. 
12 THE MAIN PART OF THE BOUNDARY LAYER_. Z3 THE_TH INNER ADJUSTMENT 
REGION CONTAINING THE CRITICAL LAYER ZA, AND Z5 IS THE WALL LAYER. 
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M^Af (y-1)/6 c BS |AJ 

(a) PLOT OF THE NEUTRAL WAVESPEED C Q AS A FUNCTION OF 
THE SCALED MACH NUMBER . Cg S IS THE VALUE OF C Q 

WHEN = 0. i.e. THE SMITH & BODONYI (1980) CASE. 
THE DOTTED LINE IS THE INVISCID LIMIT FOR INCREASING 
M^ , SEE SECTION (2.19). 



(b) THE SOLID LINE IS THE NEUTRAL PHASE SPEED c FROM 
THE GENERALIZED INFLEXION POINT CRITERION AND THE 
DASHED LINE THE ASYMPTOTE (2,19) REWRITTEN FOR (1) 
MACH NUMBERS. SEE SECTION 2. 

FIGURE 5. 



(a) THE COMPUTED PHASE SHIFT it AS A FUNCTION OF A c 
FOR VARIOUS VALUES OF THE PARAMETER J IN EQUATION 
(3.10). THE SOLID CURVE IS THAT FOR J = 0. i.e. THE 
HABERMAN (1972) CASE. 



(b) THE VORTICITY JUMP H + -H~ AS A FUNCTION OF A c FOR 
THE SAME VALUES OF J AS IN FIGURE 9 . 


FIGURE 9. 
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(a-c) A SKETCH OF THE POSSIBLE NEUTRAL CURVES FOR THE 
NONLINEAR THEORY. A„ IS THE NEUTRAL AMPLITUDE, 0 Q 

IS THE NEUTRAL FREQUENCY. L,R,BSG DENOTE THE LINEAR 
RAYLEIGH, AND BODONYI , SMITH & GAJJAR (1983) LIMITS 
AND SCALINGS RESPECTIVELY, (a) J<J C , ) J c < J AND 
^crit < A c, (c) J c < J AND 0 < A c < A c rit. 

FIGURE 5. 
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FIGURE 6. - THE PHASE SHIFT ® 1 CA c > AND THE VORTICITY FIGURE 7. - A COMPARISON OF THE PHASE SHIFT FROM THE FULL 

JUMP (fl + -H") FOR THE CASE OF |J| >>1. SEE SECTION 4. NUMERICAL SOLUTIONS Of (3.10) WITH THE ASYMPTOTE FOR 

LARGE J. AN, BN DENOTE THE NUMERICAL SOLUTIONS FOR 
\ = 1,0.5 RESPECTIVELY, AND AA, BA THE CORRESPONDING 
ASYMPTOTES FOR LARGE J. 



FIGURE 8. - A COMPARISON OF THE CALCULATED PHASE SHIFT 
WITH THE ASYMPTOTE (4.8) FOR SMALL VALUES OF A c , FOR 
DIFFERENT VALUES OF J. THE DASHED CURVES ARE THE 
ASYMPTOTES AND THE SYMBOLS ARE THE PHASE SHIFTS FROM 
THE NUMERICAL COMPUTATIONS. 
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layers and Blasius flow are considered with and without heat transfer, and the neutral eigenrelations incorporating 
compressibility effects are obtained explicitly. The compressible nonlinear viscous critical layer equations are 
derived and solved numerically and the results indicate some solutions with positive phase shift across the critical 
layer. Various lihiiting cases are investigated including the case of much larger disturbance amplitudes and this 
indicates the structure for the strongly nonlinear critical layer of the Benney-Bergeron (1969) type. Finally we 
show also how a match with the inviscid neutral inflexional modes arising from the generalized inflexion point 
criterion, is achieved. 
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